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We study the stability of topologically protected zero-energy flat bands at the surface of nodal 
noncentrosymmetric superconductors, accounting for the alteration of the gap near the surface. 
Within a selfconsistent mean-field theory, we show that the flat bands survive in a broad temperature 
range below the bulk transition temperature. There is a second transition at a lower temperature, 
however, below which the system spontaneously breaks time-reversal symmetry. The surface bands 
are shifted away from zero energy and become weakly dispersive. Simultaneously, a spin polarization 
and an equilibrium charge current develop in the surface region. 
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Introduction. The topological properties of gapless 
electronic systems have recently attracted much atten¬ 
tion m- An important example are time-reversal- 
symmetric noncentrosymmetric superconductors (NCSs) 
rnsEHni, which are characterized by strong antisym¬ 
metric spin-orbit coupling (SOC) and a parity-mixed 
pairing state [12] ■ Many NCSs display evidence of gaps 
with line nodes ffsng. This is exciting, as the line nodes 
of NCSs with dominant triplet pairing are topologically 
nontrivial defects in momentum space m- Zero-energy 
flat bands of Majorana fermions are predicted to appear 
within the projections of these nodal lines onto the sur¬ 
face Brillouin zone (BZ). Such flat bands have clear ex¬ 
perimental signatures such as sharp zero-bias peaks in 
tunneling spectra mm, equilibrium currents parallel to 
the interface between the NCS and a ferromagnet mm, 
and characteristic quasiparticle interference patterns [2T| . 

The topological properties of NCSs and consequently 
the protection of the surface states are controlled by 
the superconducting gaps, which arise from interactions. 
Properly accounting for these interactions may qualita¬ 
tively alter the surface physics. For example, a surface 
tends to suppress some gap components and enhance oth¬ 
ers [92H951 . This may change the conclusions of the afore¬ 
mentioned studies [Ml [MU, which imposed unrealis¬ 
tic uniform gaps. Flat bands with their high density of 
states are particularly prone to instabilities. Indeed, the 
zero-energy flat bands at the (110) surface of d -wave su¬ 
perconductors with time-reversal symmetry (TRS) [29] 
are predicted to be unstable towards a time-reversal- 
symmetry-breaking (TRSB) state [22H2S1 f30] . This has 
been supported by some tunneling and transport experi¬ 
ments [3Tl - f33] but was not seen in others [34U38] . d -wave 
superconductors are however qualitatively different from 
NCSs in that the zero-energy flat bands are degenerate 
in the first case but nondegenerate in the second. 
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In this paper, we study the stability of the surface zero- 
energy flat bands of nodal NCSs by performing selfconsis¬ 
tent mean-held (MF) calculations in real space for a slab 
of finite thickness. For concreteness, we consider a model 
with point group C4,,, which is realized for CePtaSi [59] , 
CeRhSia [40j, and CelrSia [41] . We show that an instabil¬ 
ity to a TRSB state can occur and study its signatures. 

Model and mean-field theory. We start from a tight- 
binding Hamiltonian for an NCS with C^ v point group, 
H = Hq + H lnt ■ The noninteracting part is 

Ho = -tY 1 C l C l ~ + C ] Ci ) 

I <*J> 

+ i A^(z x e ZJ ) ■ (c\ * Cj - cj | c*) , (1) 

<*i) 

with the chemical potential /i, the nearest-neighbor hop¬ 
ping amplitude t, and the Rashba SOC strength A. The 
SOC term breaks inversion symmetry. The annihilation 
operator Cj = (Oq, Cj-q) T is a two-component spinor, <x 
is the vector of Pauli matrices, and e,,j is the unit vector 
pointing from site j to site i of a simple cubic lattice. At¬ 
tractive interactions at the same site and between nearest 
neighbors in the xy plane are described by 

Hint = ~U S C rt C ]± C ji C jt ~ U t C \cr C \u' C jer' c i<7- 

j (ij)J Lz a a' 

( 2 ) 

The interaction is decoupled in the pairing channel. We 
define the singlet and triplet order parameters A* = 
(U s / 2) ( cJia v Cj ) and Ab = iU t ( cjicr v crci ), respectively, 
where the site indices i, j in Ab are restricted to nearest- 
neighbor sites in the xy plane. The triplet vector order 
parameter is taken to be parallel to the effective SOC 
field, Ab = Ab z x e„. This choice avoids the triplet- 
pair-breaking effect of the SOC, and is therefore energe¬ 
tically favorable in the bulk [22] . 

We first consider the MF solution for an extended 
system, assuming spatially uniform gaps A® = A s and 
Ab = A t . Details of the calculation are given in Sec. I of 
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FIG. 1. (Color online) (a) Projection onto the (101) plane 
of the positive-helicity Fermi surface (thin solid line), the 
negative-helicity Fermi surface (dashed line), and the super¬ 
conducting gap nodes on the former (heavy solid lines), for the 
bulk NCS. The gray areas denote the zero-energy flat bands 
predicted to exist at (101) surfaces under the assumption of 
uniform gaps PE]. The plot is restricted to momenta in the 
(101) surface BZ, where k m = ( k x — k z )/\/2. The parame¬ 
ters are t = 1, A = 1.5, fi = —3, U a = 5.0, Ut = 5.4, and 
T = 0.0025. (b) Mean-field gaps A s (solid black) and A t 
(dashed red) as functions of temperature T. 


the Supplemental Material [35]. We find that the singlet 
and triplet gaps have the same phase, which can be set 
to zero, so that TRS is preserved. SOC splits the bands 
and thus also the Fermi surface according to the helicity 
of states [2] ■ Since the triplet order parameter is parallel 
to the SOC, pairing only occurs between states with the 
same helicity. We determine interaction strengths U s , Ut 
that lead to flat zero-energy surface bands under the as¬ 
sumption of uniform gaps. The resulting surface states 
have been studied in detail in Refs. gum 121- This is 
realized for the parameters t = 1 (hence, t is our unit of 
energy), A = 1.5, p = —3, U s = 5.0, U t = 5.4 at the tem¬ 
perature T = 0.0025 (setting ks = 1), giving bulk MF 
gaps A s = 0.704 and A t = 1.006. We consequently find 
a gap with line nodes on the (smaller) positive-helicity 
Fermi surface, but a full gap on the (larger) negative- 
helicity Fermi surface [15]. Figurejlja) shows the projec¬ 
tion of the two Fermi surfaces and the nodal lines onto 
the (101) plane. The topological argument from Refs. 
mm predicts that a ( 101 ) surface hosts flat zero-energy 
bands within the region bounded by the projected nodal 
lines. In addition, there is an arc of zero-energy states 
connecting the two regions with flat bands mmm- Fig¬ 
ure [ljb) shows the bulk gaps A s and A t as functions of 
temperature. 

We next turn to the MF solution for a slab of thickness 
W with (101) surfaces. We introduce new coordinates 
x = m + (l + l mod 2)/2 and z = — m + (l — l mod 2)/2, 
where m is parallel to the surfaces and l = 0,..., W — 1 
is orthogonal to them. The geometry of one surface and 
our coordinate system are depicted in the inset of Fig. [2] 
Since translational symmetry in the normal direction is 
broken, the gaps depend on l. We define 
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FIG. 2 . (Color online) Selfconsistent gaps Af, Af + 1 y 2 , A? for 
a slab of thickness W = 300 and parameters as in Fig. [l] (a) 
and (b) show the real and imaginary parts, respectively. The 
lines denote the bulk gaps A s (solid black) and A t (dashed 
red). Inset: Sketch of the bottom (l = 0) surface of a (101) 
slab, showing the new coordinates l and m. The y axis points 
into the plane of the drawing. 


iUt (cjia v crci ) 


A?+i /2 z x for x bonds, 
A f z x e.y for y bonds, 


(4) 


where the subscript l denotes the (identical) l coordinate 
of sites i and j, while l + 1/2 in A) c + 1 y 2 is the mean 
of the l coordinates of sites i and j. We Fourier trans¬ 
form in the directions parallel to the slab, introducing 
the two-dimensional momentum vector k = (k m ,k y ) in 
the surface BZ, defined by —n < k y < n and — 7 r /\/2 < 
k m = (k x - k z )/\/2 < n/V2. The MF calculations are 
performed for a slab of thickness W = 300, using the 
same parameters as for the bulk calculation. Further 
details are presented in Sec. II of the Supplemental Ma¬ 
terial [43] . 

Spontaneous breaking of TRS. Our central results are 
summarized in Figs. [2] and [3] at sufficiently low tem¬ 
peratures, the singlet and triplet gaps develop imaginary 
components close to the surface, spontaneously break¬ 
ing TRS. This solution is degenerate with a state with 
complex-conjugated gaps. In the limit W —> oo, the two 
surfaces are decoupled and there are hence four degener¬ 
ate TRSB solutions, differing in the signs of the imagi¬ 
nary parts of the gaps close to the surfaces. 

The spatial variation of the gaps near the surface in the 
TRSB phase is shown in Fig. [2] While both the singlet 
and triplet gaps develop imaginary components near the 
surface, the real parts of the singlet and triplet gaps are 
enhanced above and suppressed below their bulk values, 
respectively. The suppression of the triplet gaps origi- 














3 



FIG. 3. (Color online) Selfconsistent gaps Af, Af +1 , 2 , A/ for 
the surface layer (l = 0, filled symbols) and at the slab center 
(.I = W/2 — 1, open symbols) as functions of temperature. The 
thickness is W = 300, the parameters are as in Fig. [I] (a) 
and (b) show the real and imaginary parts, respectively. The 
imaginary parts for l = W/2—1 would be indistinguishable 
from zero and are omitted. The lines in panel (a) denote the 
bulk gaps A s (solid black) and At (dashed red) from Fig. Se¬ 
nates from the pair-breaking effect of the surface, which 
in turn enhances the singlet gap to compensate for the 
lost condensation energy. The reversal of the suppres¬ 
sion of the triplet gaps in the outermost layer can be 
understood similarly: since one of the triplet amplitudes 
is missing at the surface, the others are enhanced. 

The gaps converge to their bulk values as we move 
away from the surface; the gaps at the center of the slab 
are within 0.01% of their bulk values. Note that the 
deviation of the imaginary parts from their bulk value 
(of zero) has a much longer range than that of the real 
parts. Indeed, close to the center of the slab, we find 
that Irn A/ oc (l — W/ 2), see Fig.[2](b). We have checked 
that the proportionality constant decreases more rapidly 
than W~ 1/2 with W so that the gradient energy vanishes 
for W —> oo. We attribute the slow spatial decay to the 
enhancement of length scales close to the bulk quantum 
phase transition to a nodeless singlet-dominated state. 
This transition can be reached by increasing U s and de¬ 
creasing XJt by only 0.067 (not shown). 

The evolution of the TRSB state with temperature is 
shown in Fig. |3j where we plot the gaps Af, A* +1 / 2 , and 
a r in the surface layer and at the slab center. Upon in¬ 
creasing the temperature, the gaps in the surface layer 
show a second-order transition, at which the imaginary 
parts vanish and TRS is restored. This occurs at a tem¬ 
perature of T s « 0.083, well below the bulk supercon¬ 
ducting transition temperature T c ss 0.942. 

Dispersion and density of states. In Fig. [4ja) we plot 
the dispersion for a cut through the surface BZ at k y = 0 
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FIG. 4. (Color online) (a) Dispersion for a cut through the 
surface BZ at k y = 0, for W = 300 and the same parameters 
as in Fig. [T] The black points refer to T = 0.0025 -C T s in 
the TRSB state, whereas the cyan (light gray) points in the 
background refer to T = 0.1 > T s with restored TRS. The 
dispersion is odd in fc m , only points for k m > 0 are shown, 
(b) Surface DOS in the l = 0 layer at the same temperatures. 
An artificial broadening of r/ = 0.01 was used. 


at temperatures below and above T s . For T > T s , the 
zero-energy flat band predicted in Refs. mm is clearly 
visible for 0.5 < k m < 1.5; the zero-energy states at 
k m /S 0.5 form an arc connecting the projections of the 
nodal rings [2j QT]. The TRSB for T < T s removes 
the topological protection of the zero-energy flat bands 
of the TRS state, which are consequently pushed away 
from zero energy, with a low-temperature energy shift 
on the order of T s . Since the shift is weakly momentum 
dependent, the band obtains a nonzero velocity. Due 
to particle-hole symmetry, the dispersion is odd in k. 
The zero-energy flat bands give a singular contribution 
to the surface density of states, which can be detected as 
a sharp zero-bias peak in the tunneling spectrum of an 
NCS-normal-metal junction jU [TUI m • The shift of the 
surface bands in the TRSB state causes a splitting of this 
peak, as shown in Fig.|4j(b). This splitting is a key experi¬ 
mental signature of TRSB. Indeed, the observed splitting 
of the zero-bias peak for tunneling into the (110) surface 
of the cuprates is important evidence for TRSB in this 
system j5T, ,52]. 

Spin polarization. Broken TRS is also manifested by a 
nonzero spin polarization near the surface, which is di¬ 
rected along the y-axis. A polarization in other directions 
is forbidden by mirror symmetry in the xz plane. Fig¬ 
ure [5^a) shows the spatial variation of the layer-resolved 
spin contributions (s^); explicit expressions for the spin 
operator s i in layer l and its thermal average are given 
in Sec. Ill of the Supplemental Material |43| . It is in¬ 
teresting to examine how states at different k contribute 
to the spin polarization: due to the strong polarization 
of the flat-band surface states in the TRS state [44j |46], 
one might expect that the spin polarization largely orig¬ 
inates from the shifted flat bands. To check this, we plot 
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FIG. 5. (Color online) (a) spin polarization ( sf) (filled circles) and current (j™) (open squares) as functions of l, for W = 300 
and the same parameters as in Fig. [I] Both quantities are given in units of their value at the surface. The layer indices of 
(j ™) are given as half integers to indicate that the current flows between two layers, see the inset of Fig. [5] (b) Momentum- 
resolved contributions to the y component of the total spin polarization of half the slab (0 < l < W/2) in the surface BZ. (c) 
Momentum-resolved contributions to the m component of the current in half the slab (0 < l < W/2). The momentum-space 
plots in (b) and (c) are restricted to a region just enclosing the projection of the positive-helicity Fermi surface. 


in Fig. [§b) the momentum-resolved contribution to the 
spin polarization of the half slab defined by 0 < l < W/2 
[4,*)) . Surprisingly, the spin polarization is not primarily 
carried by the shifted flat bands but rather by bulk and 
perhaps dispersing surface states HE] from the region 
between the projected nodal rings. 

Equilibrium currents. Furthermore, the absence of TRS 
permits a nonzero equilibrium surface current |22ll24[ l26]. 
Indeed, we expect such a current since the surface bands 
become dispersive and the dispersion is odd in k m \ a sim¬ 
ilar modification of the electronic structure at an inter¬ 
face with a ferromagnet does result in a surface current 
nanD]. Explicit expressions for the current operator j; 
in layer l and its thermal average are given in Sec. IV 
of the Supplemental Material [43j. Although charge is 
not conserved in the superconducting MF state, one can 
account for the pairing potentials by adding so-called 
source terms to the continuity equation 0Z1- For self- 
consistently calculated gaps, however, the thermal aver¬ 
age of the source terms vanishes, and charge conservation 
is retained m- This implies that the current perpendic¬ 
ular to the slab’s surface, i.e., in the l direction, must van¬ 
ish. Mirror symmetry in the xz plane forbids a current 
along the y -axis [43] , leaving only the current along the m 
direction, defined as (j™ 1/2 ) = ((j 1 + 1 / 2 ) ~ 0 'f+i/ 2 »/^ 2 - 
O'J+ 1 / 2 ) indeed nonzero in the TRSB state: In Fig. 
[5ja) we plot the current as a function of the layer index 
l , which shows that it is bound to the surface with spa¬ 
tial profile similar to the spin polarization. In contrast 
to the spin polarization, the main contribution to the 
current stems from surface states within the projected 
nodal rings, as shown by the momentum-resolved cur¬ 
rent in a half slab plotted in Fig. [5jc). We have also 
studied the contributions to the vanishing components 
(j\) and Of), shown in the Supplemental Material [43] , 
Interestingly, (jj) cancels only in the sum over the full 
surface BZ, showing that bulk states must be included to 


satisfy charge conservation. Note that the sign of both 
the spin polarization and the current is reversed for the 
degenerate solution with complex-conjugated gaps. 

The coupling to the electromagnetic field, which is not 
included here, leads to additional screening currents that 
exactly balance the spontaneous surface current in the 
limit W —> 00 . However, these currents build up on the 
length scale of the magnetic penetration depth A, which 
in typical NCSs is much larger than the decay length of 
the surface current, on the order of the coherence length 
£ m■ In samples with thickness smaller than the pen¬ 
etration depth but larger than the coherence length, it 
should thus be possible to detect the surface current. 

Summary and conclusions. We have studied the sta¬ 
bility of zero-energy flat bands at the surface of an NCS 
within selfconsistent MF theory. We find that the flat 
bands are indeed recovered by the selfconsistent calcu¬ 
lation within a broad temperature range below the bulk 
transition temperature T c . TRS is spontaneously bro¬ 
ken at a much lower temperature T s , which is signaled 
by a nonuniform phase of the gaps. This destroys the 
topological protection for the flat bands, shifting them 
away from zero energy and giving them finite velocity. 
Figure [3] shows that at low temperatures the flat bands 
are displaced by an energy on the order of T s , which is 
significantly smaller than the bulk gaps of order T c . The 
free energy gain due to the shift of the flat bands is likely 
a major driver of the TRSB state, and ultimately limits 
T s as the free energy gain from the shift is reduced by 
the broadening of the Fermi function. 

The TRSB state leads to clear experimental signa¬ 
tures: a splitting of the zero-bias peak in the tunnel¬ 
ing spectrum, a nonvanishing spin polarization at the 
surface, and a nonvanishing equilibrium charge current 
parallel to the surface. The latter two effects show that 
the TRSB state found here is qualitatively different from 
that predicted for the (110) surface of cuprate supercon- 
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ductors (22HISI [30] • 
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I. MEAN-FIELD THEORY FOR THE BULK 


In this section we sketch the MF theory for the bulk NCS. We assume spatially uniform pairing potentials A* = A s 
and AT = A t . Using this ansatz to decouple the interaction Hamiltonian H- mt , we obtain the Bogoliubov-de Gennes 
(BdG) Hamiltonian jSl] 

ffMF = ^^(k)$ k + Y^+iv|i, (SI) 

k 

with the number of sites, N, and the block matrix 

n(k) - ( Mk) A(k) ^ rsoi 

V At (k) -h T (-k) J (S2) 

written in terms of h( k) = £ k a-° — Al k • <x, A(k) = (A s cr° -(- A t l k • cr) ia v , £ k = —2 1 (cos k x + cos k y + cos k z ) — /r, 
l k = x sin ky — y sinfc x , and the Narnbu spinor <f> k = (c k f-, C 14 , c_ k cy k ^) T . Here, <r° is the 2x2 identity matrix. 
The dispersion U k|y , v = 1,..., 4 is obtained by diagonalizing "H(k). A s and A t are then obtained by minimizing the 
free energy 
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where the momentum sum is over half the BZ, k m > 0. This restriction of the sum makes use of particle-hole 
symmetry, which relates the Hamiltonian in Eq. (S2| at k and —k by [ST] Uc'H T {— k)tlf, = — 'H(k) with the unitary 
matrix Uc = a x ® cr 0 . 


II. MEAN-FIELD THEORY FOR THE SLAB 


We now set up the MF Hamiltonian for the (101) slab and describe the determination of the gap parameters Af, 
A?+i/ 2 , and Af in the MF approximation. After Fourier transformation in the directions parallel to the surfaces, the 
MF Hamiltonian reads 


. w-i . w -2 

H mf = X J2 $ L^«( k ) $ kz + ~ J2 J2 $ L+I%d( k ) $ : 
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(S4) 


A r w-l at W — 2 at W-l 

s 1=0 1 1=0 1 1=0 

where A) is the number of unit cells of the slab and <l > k ; = (c k ^, c k 4 , c_ k l c_ k ; , ) T is the partially Fourier- 

transformed Nambu spinor. The sums over l containing $ k ;±1 are restricted in such a way that Z± 1 £ {0,..., W — 1}. 
The coefficient matrices appearing in F/mf are 
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We next construct the 4 W x 4 W block matrix 


n k) = 


( %oo(k) 
’Hio(k) 
0 

V ; 


ftoi(k) 0 

'Hn(k) '^i 2 (k) 

w 21 (k) n 22 { k) 


(S7) 


and denote its eigenvalues by E k„, z/ = 1,... ,4W and the corresponding eigenvectors by | kzx). The MF Hamiltonian 
satisfies particle-hole symmetry [ST], WcH r (- k)k/^ = —Hik) with the unitary matrix Uc = 1 w <8> (? x <8 cr°, where 
1 w is the W x W identity matrix. This symmetry again allows to restrict the momentum sums to half the BZ. The 
free energy can then be written as 


Fmf = —ksT ^ In f 2 cosh 
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where the momentum sum is restricted to half the BZ, k m > 0. Minimization of Fmf gives the gaps Af, A^ +1 , 2 , 
and Af. The derivatives of Fmf with respect to the complex conjugate gaps can be calculated with the help of the 
Hellmann-Feynman theorem, for example 


9F M f 

dAf* 
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The momentum sums are performed on a 50 x 50 mesh, referring to the full surface BZ. Quadrupling the number of 
points in the mesh to 100 x 100 leads to changes in the MF gaps on the order of only 0.1%. 

Solving the resulting MF equations by iteration turns out to be prohibitively slow for the required W, essentially 
because the minimum of Fmf is very shallow in some directions in the high-dimensional space of gap parameters. On 
the other hand, numerical minimization making use of the explicitly known gradient is reasonably efficient. We use 
the Broyden-Fletcher-Goldfarb-Shanno method implemented in Numerical Recipes [S2]. It requires an initial guess 
for the inverse Hessian. When we scan over ranges of temperatures, we use not only the converged values of the gaps 
but also the best approximate inverse Hessian from one step as starting values for the next, which significantly speeds 
up the convergence. We assume that the method has converged when no real or imaginary part of any gap parameter 
changes by more than (double) machine precision in the last step. 

For certain parameter values, we find nonvanishing gradients of the phases of the order parameters in the l direction, 
normal to the surfaces. Specifically, we find four metastable solutions, which are mapped onto each other by inverting 
the phase gradients at one or both surfaces. In the limit W —> oo, the four solutions are degenerate. For finite W, 
they split into two degenerate pairs with phase gradients that are even and odd, respectively, under reflection at the 
center of the slab. We here choose a solution with even phase gradients since then the selfconsistent solution ensures 
that the phases of Af, Af +1 ^ 2 , and Af become equal at the center of the slab; equal phases of all gaps at the center 
are expected since the bulk MF solution has equal phases. By a global phase change we can then make the phase of 
all gaps zero at the center. The phases and imaginary parts of the gaps are then odd under reflection at the center. 
Finally, of the two remaining solutions differing in the sign of the imaginary parts of the gaps, we select the solution 
with Im Aq > 0 for definiteness. The other solution leads to inverted spin polarizations and currents. 


III. SPIN POLARIZATION 

Here, we present expressions for the spin polarization. The operator of the spin per site, averaged over the directions 
parallel to the surfaces, is 

Si = iE c ^ c k!. (S10) 

II k 

Using particle-hole symmetry, the thermal spin average can be written as 

( s i ) = zTKU El tanh^=p (kv\P u <g> f ^ _^. T ) IM, (Sll) 

#1 ku ' ' 

where Pw is a W x W matrix with the components (Pw)nn' = SinSi'n'- We also consider the momentum-dependent 
contributions to the spin polarization of the half slab defined by 0 < l < W/2. These contributions are obtained by 
summing (s;) over 1 = 0,..., W/2 — 1 and removing the factor 1 /Aj and the momentum sum. 
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IV. EQUILIBRIUM CURRENT 


The second observable of interest is the current. The operators j“- denote the electron-number current from site j 
to its nearest neighbor i in the a = x,y, z direction. They can be read off from Hq in Eq. (1) in the main text, 
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The interaction term H mt conserves charge locally and therefore does not contribute to the current operator. After the 
MF decoupling, the anomalous terms do not conserve charge—they describe creation or annihilation of two electrons 
either at the same site or at neighboring sites. Such processes do not lead to currents but do introduce a source term, 
which is discussed in the main text. We average the current over layers parallel to the surface, taking into account 
that jfj and jf connect adjacent layers, whereas jA describes a current within a single layer. We then obtain the 
thermal averages, again using particle-hole symmetry, 
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where Oj+1/2) denotes currents connecting layers l and 1 + 1. The components with respect to the slab coordinates 
are 


.1 . _ O7+1/2) d” O7+1/2) , , m . _ O/+1/2) Of+1/2) 


O/+1/2) — 


02 


O/+1/2) = 


02 


(S18) 


We note that (jf) vanishes for any choice of gap parameters for our model, even non-selfconsistent ones. This is based 
on mirror symmetry in the xz plane. The current in the y direction changes sign under this symmetry operation and 
thus vanishes. 

The momentum-dependent contributions to the current in the half slab 0 < l < W/2 are obtained by summing 
(j/) over l = 0,..., Wj2 — 1 and removing the factor 1 /A) and the momentum sum. The momentum-resolved m 
component, which sums to a nonzero current, is shown in Fig. 5(c) in the main text. We present the momentum- 
resolved y and l components in Fig. |S1| The y components chancel by symmetry, as noted above. The cancelation 
of the l components, which is required by charge conservation, is only ensured for selfconsistent gaps |S3I . Large 
positive contributions from bulk states within the projected (small) positive-helicity Fermi surface are canceled by 
small negative contributions from the flat bands and from bulk states within the projected (large) negative-helicity 
Fermi surface. This shows that the bulk states must be included to satisfy charge conservation. 
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FIG. SI. Momentum-resolved contributions to (a) the y component and (b) the l component of the current in half the slab 
(0 <l< W/2). 










